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I Abstract. In this paper we develop an artificial initial boundary value problem 

04 ' for the high-order heat equation in a bounded domain Q. It is found an unique 

I classical solution of this problem in an explicit form and shown that the solution of 

D ■ the artificial initial boundary value problem is equal to the solution of the infinite 

pL^ I problem (Cauchy problem) in il. 
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1. Introduction 



When computing the solution of a partial differential equation in an unbounded 
domain, one often introduces artificial boundaries. In order to limit the computa- 
tional cost, these boundaries must be chosen not too far from the domain of interest. 
Therefore, the boundary conditions must be good approximations to the so-called 
I transparent boundary condition (i.e., such that the solution of the problem in the 

■ bounded domain is equal to the solution in the original domain). 

One of the numerical methods for the solution of problems in unbounded domains 
. is the Dirichlet-to-Neumann (DtN) Finite Elements Method. Its name comes from 

the fact that it involves the nonlocal Dirichlet-to-Neumann (DtN) map on an artifi- 
cial boundary which encloses the computational domain. Originally the method has 
been developed for the solution of linear elliptic problems, such as wave scattering 
problems governed by Helmholtz equation or by the equations of time-harmonic 
elasticity. Recently, the method has been extended in a number of directions, and 
?— ( . further analyzed and improved, by the D. Givoli [1] and others. In [I]- [5] it can be 

find more references to various numerical methods associated with solving similar 
problems. 

This question is of crucial interest in such different areas as geophysics, plasma 
physics, fiuid dynamics (see [B] - [TO])- 

The need for practical transparent (artificial) boundary conditions combining ef- 
ficiency and simplicity is evident. Such conditions must satisfy several criteria: (i) 
The resulting initial boundary value problem should be unique and stable; (ii) the 
solution to the initial boundary value problem should coincide or closely approxi- 
mate the solution of the infinite problem on the definitional domain of the boundary 
value problem; and (iii) the conditions must allow for an analytical solution or an 
efficient numerical implementation. 

In this paper we consider artificial boundary conditions for the high-order Cauchy 
problem for the heat equation. The conditions satisfy the above-mentioned criteria 
(i) , (ii) and (iii) . Similar results were taken for the Laplace equation in [11] and for 
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high-order Laplace equation (polyharnionic equations) in |12| I13j . And a transfer 
of the Sommerfeld radiation condition to a boundary of bounded domains for the 
Helmhohz operator is studied in work |14) . 

The transparent boundary condition is usually an integral relation in time and 
space between u and its normal derivative on the boundary, which makes it imprac- 
tical from a numerical point of view. Alternatively, the requirement for boundary 
conditions can be avoided when the solution of a partial differential equation is 
approximated in the form of convolution in space and time with the fundamental 
solution. An efficient approximation of this type for the heat equation is proposed 
in [15]. 



2. Main results 

In what follows, we shall need in definitions of the Holder spaces (see, e.g., [16] . 
pp. 33-34 and 117-118). 

On a cylindrical domain Q, = Q x (0,T), where Q C W\n G N is a simply- 
connected bounded domain with a sufficiently smooth boundary dQ, we consider 
the following heat potential 

(1) U{x,t) = / / EmA^ - ~ T)f{^,T)d^dT, 



Jo JQ 

where em,n{x,t) = exp(— ■^^),m £ N is a fundamental solution of the 

Cauchy problem for the high-order heat equation, i.e. this solves the following 
high-order heat equation 

and its adjoint 

iOlrre^A^ -^,t-T) = {-^~ A^re^Ax - e, t - r) = 

for alH > r and £,,x E i?". 

Indeed, it's not difficult to check that for m > 1 

<>x,t£7nA^ - ^, i - t) = em-i,n(a; -^,t-T) 

and 

<>^,r^mA^ - ^ - t) = Eni^iA^ -^t- t) ■ 

From [17, chapter 1] for x <E Q the condition 

lim / (0+J"-i£^_„(a;-^,t-r)/(e,r)dC = /(x,i) 

is valid. ^ _^ 

It is easy to see that if f{x,t) G C'2'/ {^) then u{x,t) € C^.™^'''™^^ (fi), where 
< 7 < 1 and 

(2) <}"'u{x,t)^f{x,t),{x,t)€n, 
with initial conditions 

d'^uixA) , „ 

(3) ' \t=o^O,t^O,m^l. 
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We note that if the function u{x,t) G (^^™+'''^™+2 j-j-j^ satisfies the conditions (|3]) 
then it will satisfy to the following conditions too 



<>*-u(a;, t) |t=o= 0, i = 0, TO - 1. 
Consider the problem ([21)-© with non-local type boundary conditions 

IA\ rki 4-\ I — ()^u{x, t) 

(4) Iu[x,t) UeaQ= 

+ E [/ / if^ -0'''-'-'u(i,r)dS^dT 

Jo JdQ 9n^ 

— ^-^{OU+'SmA^ - - T)dS^dT] = 0, 



JdQ 



fc = 0, m - 1, (x, t) edQ X (0, T), 



where jf— denotes the exterior normal derivative on the boundary dQ 



Theorem 1. The heat potential for the high-order heat equation (0) is an unique 
classical solution of the non-local type initial boundary value problem 

Proof. The following equalities are valid 

(❖^+J"e™,„(x-e,i-T) = (❖+J"-^£™-i,„(x-C,i-T) = ... = Ol,e,,n{^-^,t-r) = 0, 
for alH > r and e i?", 

Iini(0+j'^e™,„(.x-e,i-^) =0, 

for fc < m — 1. Therefore 



lim / (O+^)''£„,„(x-C,a)u(C,t-aK = 



for fc < m — 1 and 

' {0+rem.n{x^^,t~r)ui^,r)d^dT^ lim / I iOt,rem.n{x-^,t-T)ui^,T)d^dT 



lim / / Ou{^,T)d^dT = 0. 
=-^0+0 io Jo 



We assume that u{x,t) G C^T^"''''''™^^ (fi). Let us denote 0" = /, where / is 
identity operator. A direct calculation shows that 

ft— a /• 

u{x,t) ^ lim / / e„^.n{x -^,t-T)f{^,T)d^dT 

= lim / f eraAx-tt-T)<^{<}"'-^u{tT))d^dT 

or 



m — 1^ 



JQ 

t-Q 



a-i-O+O 



JQ 



JdQ dn^ 



/ / em.nix ~ t,,t - T)dS^dT\ = ... 

Jo JdQ dn^ 

= lim [/ {0+)"'-h^.n{x-^,t-T)u{tT)\l" d^ 
* " / iOtrrem.nix ~ - T)uitT)didT] 

Jq 

+ $:[// ^""^"^ "'2 -0'''-'-M^,r)dS,dr 

~5 Jo JdQ on^ 



d^^ 



-i<>try^rn.nix-(,t~T)dS^dT] 

/O JdQ 



/o Jag 



for any {x, t) £ fi. 
Here we get 



Jo JdQ 

^^-^^ ^(0+ re„,„(a;-C,i-T)d^5dT] =0,V(x,i) ef^. 

Applying properties of the double layer and single layer potentials (see |17] or 
[in]) to © as {x,t) dQ X {0,T), we obtain 

2 ^ Jo -/OQ 



10 JdQ 



{<>l,ye„^,n{^-C.t-r)dS^dT]=0,y{x,t) edQ x (0,r). 



Applying the differential expression ^ , fc = 1, m — 1 to u{x, t), we get 



t-Q 



Va^O+Ojn 



lim 







= 0./ [J_lZoJ„ j^e„,.iA^-i,t-T)f{^,T)d^dT 

nt — Ct p 

^^"^.n / ^'m-k,n{^ - ^ - T)d^dT 



"^0+0 JO JQ 

and after similar calculations as above, we get the following boundary conditions 
for each k 

= ^ ^+ 2^ / / 0" ^u{tT)dS^dT 

2 -'o Jag 

dn^ ^' 



V(a:, t) edQ X (0, T), /c = 0, TO - 1. 

Thus, the heat potential ([T]) satisfies the boundary conditions (H]). 

Conversely, if a function ui{x,t) G cf 7'*'''^'"'''^ ^ (^) satisfies equation ([5]), the 
initial condition ([3]) and boundary conditions Q then ui{x,t) = u{x,t), where 
u{x, t) is the heat potential ([T]). If this is not so, then a function v{x, t) ~ ui(x, t) — 
u{x,t) satisfies the homogeneous equation 

(6) <)"'v{x,t)=0,{x,t) eil, 
with the initial conditions 

(7) <>'"v{x,t) \t=o= 0,x e n,k^ 0,TO- 1, 
and the boundary conditions 

(8) -^^'(2^,*) = 0,V(x,t) edQx (0,r),fc = 0,TO- 1. 
Where we use the following notations 

(9) I^ix,t)^"^f:\f [ ^^^^^-^'^''7^""^''""^ 0— ^t^(^,r)dg,dr 

^0 JdQ dn^ 

JdQ on^ ^' 



V(x,i) e 17,fc = 0,m - 1. 
On the other hand, by using ([B]) and ([7]), we obtain following equalities 



= / / e,n,nix -^,t-T) ■ Od^dr 
Jo JQ 







= 0'=w(a;, + I^{x, i), V(x, G 17, fc = 0, to - 1. 
Applying properties of the double layer potential (see [17], [18] or [19]) to dU as 
(x, t) ^ dQ X (0, T), we have 

\dQx(o,T)^ ~Iv{x,t) laQx(o,T)= 0,fc = 0,m - 1. 
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I.e. the problem (l6])-([8]) is equivalent to 

(10) =0,(x,t) e 17, 
with the initial conditions 

(11) <>''vix, t) |t=o= 0, X e Q, fc = 0, m - 1, 
and the boundary conditions 

(12) <?''v{x,t) |aQ^(o^r)=0,fc = 0,m-l. 

From ([T0| -(fT2 |) with k = m ~ \ using the Maximum Principle we get 

<?'"'-^v{x,t) = 0,(a;,t) e VI. 

And repeating the similar calculations we have 

C>''w(a;, t) = 0, V(x, t) eVl,i = 0,m- 1, 

i.e., w(a;, i) = Ui(a;, t) — M(a;, t) and Ui{x, t) = u{x, t). 

This completes the proof of Theorem [T] □ 



Remark 1. The heat kernel, i.e. the fundamental solution of the high-order heat 
equation em,n{x,t) is the Green function for the non-local type boundary value prob- 
lem 

Let us consider the heat equation 

(13) <>u{x,t)^ f{x,t),{x,t)en, 
with initial condition 

(14) u{x,Q) = Q 

and with inhomogeneous non-local type boundary condition 
u{x,t) /■* f dei,n{x - £.,t - t) 

(15) — + / / U(e, T)dS^dT 

^ Jo JdQ on^ 

f ^^(^ 

- — ^ — ei^nix - £,,t - T)dS^dT = ip{x,t), 

Jo JdQ Cl^i 

{x,t) edQx (o,r), 

where ip{x,t) € cf ^'^'^"^^ (9Q x [0,r]), <p(x,0) ^ 0,x e dQ and denotes the 
exterior normal derivative on the boundary dQ. 



Remark 2. In case ip{x,t) = 0, from the theorem[^ follows that an unique solution 
of the problem !il3\} - ![T^} is given by potential 

(16) U{x,t)^ f f ei,n{x^^,t-T)f{tT)d^dT. 

Jo Jq 
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Lemma 1. The function 

17) U{x,t)= / ip{^,T)d^dT 

Jo JdQ dn^ 
is an unique solution of the following problem 

(18) <}u{x,t) =0,{x,t) en, 
with initial condition 

(19) u(x,0)=0 

and with inhomogeneous non-local type boundary condition 

u{x,t) /■* f dei,n{x - £.,t - t) 

20) — + / / j: u{£„T)dS^dT 

2 Jo JdQ dn^ 

-II ^"i^' e,n{x -^,t- r)dS^dT = ^{x, t), 
Jo JdQ 

{x,t) edQx (0,T), 

where denotes the exterior normal derivative on the boundary dQ and G{x, ^, t, r) 
is the Green function of the problem JJ3)) -( |T^ with the Dirichlet boundary condi- 
tion. 

Proof. 

= I I ei,n{x - £„t - T)Qd£,dT = f f £i,„(x-C,t-T)Ou(C,r)d^dT = 

Jo JQ JO JQ 

^u(x,t)+f f <}+e,,n(x~^,t~T)u{^,T)d^dT-ff ^^l^^-M^^i^ r)dSidT 
Jo Jq Jo JdQ on^ 

eUx-U-r)^^$^dS^dr^u{x,t)- f [ ^^lA^-i2l^u{^^r)dS^dr 
lo JdQ cln^ Jo JdQ on^ 

+ f I e,,Jx -^,t- r)^^^^dS^dT 
Jo JdQ 

for all {x,t) e g X (0,r). Hence, 

u{x,t) = / / u{^,T)dS^dT~ / ei^n[x~^,t-T) — dS^dr 

Jo JdQ dn^ Jo JdQ dn^ 

for all {x,t) e g X (0,r). And 

f ^\ u{x,t) /■* f dei.n{x - £_,t - t) 

u{x,t) = h/ / ^ u{^,T)dS^dT 

2 Jo JdQ 

~ f I e,A^ -U- T)^^$^dS^dT 
Jo JdQ 



as {x,t) dQ X (0,r). Since, 



f +\ u{x,t) f f dei.n{x - i,t - t) 

u{x,t)^ h/ / ^ u{^,T)dS^dT 

2 Jo JdQ dn^ 
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Jo JdQ dn^ 
for all e dQ x (0,T). So, problem ([T8])-(l20l) equivalent to the problem for 

homogeneous heat equation (ITSl) with initial condition (|T9l) and with the following 
Dirichlet boundary condition 

t) = ip{x, t) 

for all (x, t) G x (0, T). The lemma is proved. □ 
From lemma [T] and remark [2] follows statement of the following theorem. 



Theorem 2. The solution of the problem lll.S\) - ![T5\) is given by the following formula 
(21) 

u{x,t)= f f eU^-^,t-T)f{^,T)d^dT+ f f ^^^^'^'^'^^ y(g,r)dgdT, 
Jo Jq Jo JdQ one, 

where G{x,S,,t,T) is the Green function of the problem 1113 \) - [T^ with the Dirichlet 
boundary condition. 
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